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REGIONS  OF  STABILITY  OF  FEL  OSCILLATORS 


B.  Levush  and  T.  M.  Antonsen,  Jr. 

Laboratory  for  Plasma  and  Fusion  Energy  Studies 
University  of  Maryland 

\  ABSTRACT 

The  stability  of  single  mode  operation  of  FEL  oscillators  is 
investigated.  We  consider  two  models  of  an  untapered  FEL  oscillator. 
The  first  model  is  called  the  klystron  model.  In  it  the  FEL  interaction 
occurs  at  two  points:  a  prebunching  point  and  an  energy  extraction  point. 
In  this  model  the  nonlinear  electron  dynamics  are  solvable  exactly 
leading  to  a  complex  delay  equation  for  the  wave  fields.  The  stability 
of  single  mode  operation  can  then  be  determined  easily  as  a  function  of 
single  pass  gain,  energy  mismatch-frequency,  and  the  difference  between 
the  group  velocity  of  the  radiation  and  the  beam  velocity.  The  second, 
more  realistic  model  has  a  distributed  interaction  region  of  finite 
length.  Stability  of  single  mode  operation  in  this  device  must  be 
determined  numerically.  The  results  of  the  models  in  low  gain  regime 
will  be  compared  and  the  parameter  regimes  where  stable-single  mode 
operation  is  possible  will  be  determined. 

I.  INTRODUCTION 

Many  applications  of  free  electron  lasers  require  a  coherent  single 

frequency  output.  However,  the  FEL  interaction  region  is  large  and  in 

oscillators  many  modes  of  the  interaction  cavity  are  linearly  unstable. 

Therefore  it  is  desirable  to  know  under  what  circumstances  nonlinear 

operation  in  a  single  frequency  mode  will  occur.  This  problem  has  been 
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considered  previously  by  a  number  of  authors.  •  *  Our  work  closely 
parallels  that  of  Refs.  (1)  and  (2).  Specifically,  we  consider  two  cases 
in  which  reduction  in  computation  is  achieved  by  considering  a  model 
system  or  by  exploiting  a  small  parameter. 

We  begin  with  the  one  pendulum  equation  describing  the  space  and 
time  evolution  of  the  particle  energy  and  phase  and  signal  field  in  an 
FEL, 

(ft  +  vz  57^  ’  *  {Y(wvz)  fy-  (v"1)  ,  (1) 
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where  f  and  5y  are  a  particle  phase  and  energy  mismatch,  and  A(z,t)  is 
the  complex  amplitude  of  the  signal  field.  The  definition  of  the  other 
quantities  appearing  in  Eqs.  (l)-(3)  follows:  oi  and  kz  are  the  wave 
frequency  and  wavenumber,  j  is  an  effective  current  density,  v^y  )  Is 
the  mean  beam  speed,  for  the  case  considered  here  ^  is  a  constant  and 
gives  the  energy  of  a  particle  exactly  resonant  with  the  carrier  signal. 
The  undulator  parameter  is  defined  by 

qA  (z) 

..  w 


Since  we  are  considering  oscillators  it  is  necessary  to  augment  the 
wave  equation  (3)  with  a  boundary  condition  describing  the  return  path  of 
the  radiation, 

A(z=0,  t)  -  RA(z*L,  t-L/v  )  ,  (4) 

O 

where  R  is  the  combined  reflection  coefficient  for  the  two  mirrors,  L  is 

2 

the  separation  between  mirrors  and  v^k^c  /u>  is  the  group  velocity  of  the 

radiation.  (Note  that  R  is  a  complex  number.) 

We  now  make  two  transformations  to  Eqs.  (l)-(3).  The  first  is 

simply  to  normalize  all  the  variables.  Specifically  we  introduce, 
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Second,  we  define  a  new  coordinate  system  moving  with  the  mean  beam  speed 

V  -  5  (6a) 

t'  -  t  -  6*5  ,  (6b) 

where  8.*v  /v  is  the  ratio  of  the  radiation  group  velocity  to  the  mean 
£  2 

beam  speed.  With  these  changes  the  governing  equations  become, 
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with  the  boundary  condition 

a(O.T')  -  Ra(  1 ,  t'  -  (1+6*))  .  (10) 

The  function  X(?)=A  (z)/A„  gives  the  profile  of  the  coupling  coefficient 

w  wu 

between  particles  and  fields.  We  note  that  the  system  of  equations  (7)- 
(10)  has  four  dimensionless  parameters:  the  normalized  beam  current 
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where  IA=mc/q» 1.7x10  Amps,  the  relative  velocities  8*,  the  reflection 
coefficient  R,  and  the  injection  momentum  p(5'=*0). 

In  the  remainder  of  this  paper  we  will  consider  two  ways  this  system 
can  be  further  simplified. 

II.  KLYSTRON  MODEL 

In  the  klystron  raodel^’^  one  assumes  a  specific  form  for  the  profile 
of  the  coupling  coefficient, 

X(C )  -  6(5')  +  5ft'  -  1)  . 

Physically,  this  models  an  FEL  with  two  short  interaction  regions 
separated  by  a  drift  space.  In  the  first  interaction  region,  5'«0, 
particles  receive  a  kick  in  energy  (p)  depending  on  their  entrance  phase 
and  the  strength  of  the  signal  field  aft*0,T').  The  particles  then  bunch 
In  phase  as  they  traverse  the  drift  space  producing  a  coherent  current 
source  which  amplifies  the  radiation  at  the  second  interaction  region. 
The  equations  of  motion  describing  this  process  can  be  solved  exactly 
yielding  an  expression  for  the  output  radiation  field  for  a  given  input 
field.  Inserting  this  expression  in  the  feedback  boundary  condition 
produces  the  following  two  time  delay  equations. 


a(O.r')  -  r{ a(0 ,  x'-2)  +  i  |  e_lp(0) 


a(0,  T'  -  ( 1+6 *))n}  , 
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with  a  Bessel  function. 

Physically,  the  first  term  on  the  right-hand  side  of  Eq.  (11) 
represents  the  decay  of  the  circulatng  radiation  which  in  our  variables 
takes  two  units  of  time  to  complete  one  circuit  of  the  resonator.  The 
second  term  describes  the  amplification  of  the  radiation  due  to  the  beam. 
The  delay  in  this  case  is  ( 1 +0 ^ ) *  the  "one"  coming  from  the  return  trip 
of  the  radiation  and  the  "B^"  coming  from  the  time  the  particles  take  to 
traverse  the  drift  region. 

For  equal  beam  and  radiation  speeds  (6^=1)  Eq.  (11)  becomes  a 
simple  return  map  in  which  each  time  slice  of  the  radiation  field  evolves 
independently  of  its  neighbor  depending  only  on  the  value  of  the 
radiation  field  two  units  of  time  in  the  past.  This  would  also  be  the 
case  if  we  had  set  B=1  in  Eq.  (9)  without  introducing  the  klystron 

A  * 

model.  The  function  h(|a|)  in  this  case  would  be  replaced  by  some  other 
transcendental  function.  However,  as  shown  by  Fiegenbaum^  the  behavior 
of  all  such  return  maps  with  quadratic  maxima  is  universal.  In 
particular,  depending  on  beam  current  the  sequence  of  field  values 
a(0,r-2n)  may  converge  to  a  single  value  (single  mode  operation),  or  a 
periodic  sequence,  or  for  sufficiently  high  current  a  chaotic  sequence. 
We  note  that  the  return  map  model  is  only  appropriate  for  beam  currents 
low  enough  to  produce  a  sequence  of  constant  values  for  a(0,T).  This  is 
because  even  if  the  sequence  converges  to  a  period  two  or  higher  orbit 
the  function  a(0,r)  must  necessarily  have  a  discontinuity.  In  which 
case,  no  matter  how  small  1-6*  is  since  it  multiplies  the  time  derivative 
of  a(0,x)  in  Eq.  (9)  it  can  not  be  neglected.  Thus,  the  two  time  delay, 
Eq.  (11),  is  useful  in  that  it  allows  for  an  assessment  of  the  effects  of 
finite  slippage,  i.e.,  1-6^0. 

The  stability  of  a  single  frequency  nonlinear  mode  is  determined  as 
follows.  We  write  the  field  in  terms  of  an  equilibrium  value  and  a 
perturbation  * 

a(x)  -  exp(-iu>0x)[aQ  +  5a(x)] 

where  ag  is  the  constant  amplitude  of  the  single  mode  with  frequency 
and  6a(x  )~exp(-iwx )  is  the  perturbation.  The  frequency  and  amplitude  in 

equilibrium  satisfy  A  .. 

ZlW-  j  lo  n  A 

1  “  Rt  e  +  y  e  n  ( |  aQ  |  )]  >  (12) 

where  0  »(ir/2)+(l+6.)u>  -p(0).  Thus,  for  a  given  current  I,  reflectivity 
0  *  0 


R,  and  Injection  momentum  p(0),  Eq.  (12)  is  a  transcendental  relation 
determining  the  field  amplitude  8q  and  frequency  a  .  The  equilibrium 
relation  is  readily  solved  in  the  low  gain  regime  where  both  I  and  v=l-R 
are  small.  One  finds  uj  -m+6a)  with 

U  U  A 

2v  A  j  * 

—  -  cos0Qn(  1  aQ  |  )  and  5uQ  -  -  j  sin0on(  |aQ|  )  . 

Thus,  the  frequency  is  essentially  that  of  a  cavity  mode  with  a  small 
shift  proportional  to  the  beam  current,  and  the  field  amplitude  is 
determined  by  the  gain  and  the  ratio  of  the  transmission  coefficient  (v ) 
to  the  current. 

Linearizing  the  delay  equation  results  in  the  following  dispersion 
relation  for  the  frequency  u»  of  the  perturbation 

{(e*21a,-l)  -  ZR[(l+g)eik-l]}  x  { (e-21u-l)  -  ZR[elk-l]} 


(13) 


+  Z2[(l+g)eik-l](eik-l)  -  0  , 
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where  k»(8  -l)ui,  g-a^/n  (aQ)3n/3a0,  and  Z-l/2(nRIe  ),  while  ZR  and  Zj 
are  the  real  and  imaginary  part  of  Z,  respectively. 

In  the  low  gain  regime  | Z ] << 1  to  lowest  order  Eq.  (13)  gives 

(e  -  1)  -  0 

and  the  frequency  of  the  perturbation  is  also  approximately  that  of  a 
cavity  mode 


w  -  irn  (n  *  0)  . 

The  growth  or  damping  of  the  mode  is  determined  in  next  order  by  the 

eqUatf°n  .  ik  ik  ik 

W 2  -  cos0qW[ ( 2+g )e  n  -  2]  +  (e  n  -  l)[(l+g)e  n  -  1]  -  0  , 

*  -2iu)  .  .  * 

where  w*(e  — 1 ) /  |  Z  |  with  k^-ntr  (0  -1).  Regions  of  stability  in  the  k 

vs.  cos0q  plane  are  shown  in  Figs.  (1)  for  different  values  of  g.  For 
the  equilibrium  the  average  change  in  momentum  Ap  is  given  in  the  low 
gain  regime  by, 

Ap  -  -  cos0 q S a0 1  n(|a0D  . 

Thus,  the  value  of  cos0q  determines  the  dependence  of  the  efficiency  on 
the  injection  momentum.  The  quantity  g  depends  only  on  aQ  and  is  related 
to  the  derivative  of  the  efficiency  with  respect  to  the  field  amplitude, 
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The  efficiency  maximizes  for  g=*-2  and  cos9q=1. 

The  plots  in  Fig.  (la)  are  periodic  in  k  with  period  2it  (only  the 
fundamental  period  is  shown).  Thus,  for  an  irrational  value  of  (0*-l), 
knmod(2ir)  will  eventually  assume  any  value.  In  the  klystron  model, 
therefore,  only  equilibria  with  efficiency  sufficiently  maximized  with 
respect  to  cos9q  (injection  momenta)  will  be  stable.  Periodicity  with 
respect  to  k  is  an  artifact  of  the  klystron  model,  and  we  may  not  expect 
this  result  to  hold  more  generally. 

A  second  more  important  restriction  is  noted  by  examining  Fig.  (lb) 
where  g<-2,  corresponding  to  a  field  amplitude  in  equilibrium  that 
exceeds  the  value  which  maximizes  efficiency.  There  it  is  seen  that  no 
single  frequency  equilibria  are  stable. 

The  preceeding  discussion  applied  in  the  low  gain  regime  where  both 

A 

the  beam  current  I  and  transmission  coefficient  v=*(l-R)  are  small.  In 

A 

this  case  the  stability  boundaries  depend  only  on  the  ratio  I/v  through 
the  field  amplitude  aQ. 

In  the  high  gain  regime  an  additional  parameter  measuring  the 

A 

current  is  introduced,  | Z | =( l/2)RIn .  Stability  boundaries  in  the 
cos0q  -  | Z| plane  are  shown  in  Fig.  (2)  for  the  high  gain  case.  The  large 
triangle  in  the  upper  left  is  the  region  of  parameter  space  allowed  by 
the  equilibrium  condition.  The  low  gain  results  correspond  to  the  region 
of  the  plane  with  small  |z|.  As  can  be  seen  in  Fig.  (2)  where  g>-2  the 
stable  operating  regime  occupies  the  small  triangle  in  the  upper  left. 
The  stable  range  of  cos0q  values  decreases  as  current  is  raised.  Above  a 
critical  current  no  single  mode  operation  is  possible.  For  g<-2  we  again 
find  no  regime  where  stable  operation  is  possible. 


III.  MULTIMODE  SIMULATIONS 

We  now  drop  the  klystron  model  and  return  to  the  more  realistic  case 
of  a  constant  coupling  coefricient  in  the  interaction  region, 


X(C) 


{1  for  0  <  C  <  1 
0  otherwise 


To  reduce  the  number  of  parameters  we  again  assume  that  the  beam  current 
and  transmission  coefficient  are  small  and  we  are  in  the  low  gain  regime. 

A 

Solution  of  the  wave  equation  in  the  limit  1+0  will  reproduce  the  vacuum 


cavity  modes.  Including  the  current  and  transmission  losses  will  cause 
the  complex  amplitude  of  each  mode  to  evolve  slowly  in  time.  To  treat 
this  problem  mathematically  we  introduce  a  two  time  scale  formulation 
writing  the  field  amplitude 

iu> nf(l-0*)5  -t  ] 

a(5  ,t  )  58  l  an(x  )e  ,  (15) 

where  is  the  fast  time  associated  with  the  time  of  flight  of  radiation 
through  the  cavity,  to^^n  is  the  frequency,  and  t  is  the  slow  time 
associated  with  the  growth  or  decay  of  the  modes. 

The  computational  savings  is  achieved  because  the  field  is  nearly 
periodic  with  period  Tq=*2.  Thus,  an  ensemble  of  particles  can  be 
launched  with  in  the  range  0<t^<2,  the  currents  computed,  and  the  wave 

fields  advanced  by  a  time  t  of  the  order  of  the  growth  time.  This 

reduction  has  been  achieved  similarly  in  other  FEL^  and  gyrotron 

simulations .^”7 

The  governing  system  of  equations  is  found  to  be 

.  (5,Tn) 
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where  'l'n=^-^n[0-6*)5^gl  with  initial  conditions  tq  €  [0,2], 

<K0)  €  [0,2tt],  and  p(0)=*p^nj  for  the  particles.  The  wave  amplitudes  are 


advanced  on  the  time  scale  t 


where 


d  3n  1  I  .1.  ..  ^n. 

rr +  T  a„  -  -  j  (7)<<e  5 
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and  the  double  average  is  given  by 


2n  .  2  dr.  1 

«...»-/  -~f  <£(...)  . 

0  0  2  0 

A  sample  run  is  shown  in  Fig.  (3)  where  the  magnitude  of  the  amplitudes 
of  the  modes  are  shown  at  several  times.  The  injection  momentum  and 
current  to  damping  ratio  were  selected  so  that  in  the  case  of  single  mode 
operation  the  n— 3  mode  would  be  the  most  efficient.  At  early  times  one 
sees  that  the  n-*l  mode  grows  fastest  which  is  consistent  with  small 


signal  theory.  As  time  goes  by  the  system  nonlinearly  evolves  toward  a 
single  mode  (the  n=0  mode).  At  the  time  of  writing  of  this  paper  the 
code  was  still  running  so  we  can  not  report  whether  a  final  single  mode 
solution  is  reached.  However,  we  hope  to  examine  parameter  space  and 
determine  the  regimes  of  single  mode  operation  in  the  future. 

IV.  CONCLUSIONS 

We  have  outlined  two  models  capable  of  describing  multimode  effects 
in  FEL  oscillators.  In  the  klystron  model  the  wave  field  satisfies  a  two 
time  delay  equation.  In  the  low  gain  regime  single  mode  operation  is 
achieved  only  at  a  field  value  less  than  that  which  maximizes  the 
efficiency  and  then  only  for  modes  for  which  the  efficiency  has  been 
nearly  maximized  with  respect  to  mismatch.  In  the  high  gain  regime  the 
band  of  stable  mismatch  frequencies  decreases  as  current  increases.  In 
the  more  realistic  distributed  interaction  model  one  must  resort  to 
numerical  simulation.  However,  the  computation  time  and  parameter  space 
can  be  greatly  reduced  in  the  low  gain  regime. 
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FIGURE  CAPTIONS 

Fig.  1.  Stability  regions  in  the  low  gain  regime  in  cos0q  versus  k  ■ 
nrr(8  -1)  plane:  (a)  g  ■  -  1.95,  (b)  g  =*  -2.01. 

Fig.  2.  Stability  regions  in  the  high  gain  regime  in  cos0q  versus  Z  =* 
l/2(IRn)  plane. 

Fig.  3.  Sample  of  the  multimode  simulation  with  21  modes, 

g*  =•  1.2,  (I/v)  -  112.74,  and  Pinj  -  3.258;  (a) 

t  ■  1.2;  (b)  t  =*  4.8;  (c)  x  -  18;  (d)  x  *  74.4. 
s  s  s  s 


